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We study the statefinder parameter in tlie five-dimensional big bounce model, and apply it to 
diflterentiate the attractor solutions of quintessence and phantom field. It is found that the evolving 
trajectories of these two attractor solutions in the statefinder parameters plane are quite different, 
and that are difTerent from the statefinder trajectories of other dark energy models. 
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■ I. INTRODUCTION 

o : 

■ Recent observations of Cosmic Microwave Background (CMB) anisotropies indicate that the universe is flat and 
^ the total energy density is very close to the critical one with ^totai — 1 [H. Meanwhile, observations of high redshift 
O ^ type la supernovaeQ reveal the speeding up expansion of our universe and the surveys of clusters of galaxies show 

\ that the density of matter is very much less than the critical density 0]. These three tests nicely complement each 

' other and indicate that the dominated component of the present universe is dark energy. The Wilkinson Microwave 

C Anisotropy Probe (WMAP) satellite experiment tells us that dark energy, dark matter, and the usual baryonic 
matter occupy about 73%, 23%, and 4% of the total energy of the universe, respectively. The accelerated expansion 
of the present universe is attributed to the dark energy whose essence is quite unusual and there is no justification for 

, assuming that it resembles known forms of matter or energy. Candidates for dark energy have been widely studied and 

0H' focus on the cosmological constant A [3, S @, S S] with w = — 1, a dynamically evolving scalar field (quintessence) 

q! @, Elllllillillillil with m; > -1 and phantom [H [H [ll [l| with u; < -1. The fine tuning problem is 

$-H ' considered as one of the most important issues for dark energy models and a good model should limit the fine tuning 

^ " as much as possible. The dynamical attractors of the cosmological system have been employed to make the late time 

. behaviors of the model insensitive to the initial conditions qf^the field and thus alleviate the fine tuning problem, 
which has been studied in many quintessence models [13, EH, l22|, [23] ■ In phantom field, this problem has been studied 

y—{ . in [23, [2^, [2^, [231 . There are many dark energy models are constructed for interpreting the cosmic acceleration 

^ ' and solving the fine tuning problem. In order to differentiate these models, a diagnostic proposal that makes use of 

CD ', parameter pair {r, s}, the so-called statefinder, was introduced by Sahni et al.jl^ and defined as follows: 
!>■ . 

\0- a _ r-l 

' Here q is the deceleration parameter. The statefinder is a "geometrical" diagnostic in the sense that it depends on 
C . the expansion factor and hence on the metric describing space-time. Since different cosmological models involving 
' dark energy exhibit qualitatively different evolution trajectories in the s — r plane, this statefinder diagnostic can 
^ , differentiate various kinds of dark energy models. For the spatially fiat LCDM cosmological model, the statefinder 
parameters correspond to a fixed point {r = 1, s = 0}. By far some models, including the cosmological constant, 
quintessence, phantom, quintoin, the Chap lyg in gas, braneworld models, holographic models, interacting and coupling 
dark energy models [H, [2^ [s^, [3l|, [H, [H, |3J, [HI, [11] , have been successfully differentiated in the standard 4D FRW 
model. 

In Kaluza-Klein theories as well as in brane world scenarios, our 4D universe is believed to be embedded in a higher- 
dimensional manifold. In this paper, we consider a class of five-dimensional cosmological model which as an alternative 
candidate to the standard 4D FRW model, has been discussed by many authors [S^, [H, [H, [H, [HI, [H, [H, [M [H, 
[11j[13iI11]- Instead of the Big Bang singularity of the standard model, this 5D cosmological model is characterized 
by a "Big Bounce", which corresponds to a finite and minimal size of the universe. Before the bounce the universe 
contracts, and after the bounce it expands. This model is 5D Ricci-flat, implying that it is empty viewed from 5D. 
However, as is known from the induced matter theory [491. 150] . 4D Einstein equations with matter could be recovered 
from 5D equations in apparent vacuum. This approach is guaranteed by Campbell's theorem that any solution of 
the Einstein equations in N-dimensions can be locally embedded in a Ricci-flat manifold of (N+l)-dimensions [5l| . 
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The purpose of this paper is to study the statefinder parameter in the five-dimensional bounce model, and apply it to 
contrast the attractor solutions of quintessence and phantom field respectively. In section II, we introduce the scalar 
field in the 5D bounce cosmological solutions and the 5D quintessence and phantom model for dark energy, and deduce 
the statefinder parameter in 5D universe. Section III and IV are to apply the statefinder parameters to analyze the 
attractor solutions of quintessence and phantom respectively in 5D universe. Section V is a short discussion. 



II. STATEFINDER PARAMETER IN THE 5D MODEL 

An exact 5D cosmological solution was given firstly by Liu and Mashhoon in fE^ and restudied by Liu and Wesson 
in [33 . This solution reads 



dS^ = B'^{t,y)df - A^it.y) 

IdA _A 
jjL dt fi 



' {d9'^ +sm^ 9d(p' 



1 — kr"^ 



A' = (f,> + t)j,» + 2,.s + 4^. 

where A = {d/dt)A, fi = and ly = iy'{t) are two arbitrary functions of t , fc is the 3-D curvature index (fc = ±1, 0) 
and iiT is a constant. This solution satisfies the 5D vacuum equation Rab — 0, with the three invariants being 

h = R = 0, I2 = R^''Rab = 0, h = R'^''^''RABCD = ^^. (3) 



So K determines the curvature of the 5D manifold. 
The 5D line element in ^ contains the 4D one, 

ds^ = gf.^dx^dx" = B^df - A^ 



dr^ 



1 — kr^ 



r^{d0'^ +sin^ 9d(l>^) 



(4) 



Using this 4D metric we can calculate the 4D Einstein tensor by ^^^G^ =^^^ R^ ^ ^'^^R/2. Its non- vanishing 
components are 

(4)^0 _ 3 (m^ + fc) 

Wg} = (4)g2 =(4) G^ = ^ + (5) 
^ AA A^ 

Generally speaking, the Einstein tensor in ^ can give a 4D effective or induced energy-momentum tensor '^^^Tj^ via 
(4)GM = i^{i)Tif with k2 = SttG. 

In previous works [13, IH, HI, SO, SH, 113, S, H H H 113, IHl , this energy- momentum tensor was assumed to be 
a perfect fluid with density p and pressure p, plus a variable cosmological term A, and the results shows that this 
assumption works well. However, there are more and more helpful candidates for dark energy seem to be the scalar 
field models such as the quintessence and phantom. So, in Ref. [27[, the authors constructed the model that the 4D 
induced energy-momentum tensor consists of two parts: 

(4)rTi rprri 1 rp(jy 

= (Pm + Pm)u^U^ - PrnQtiv, (6) 

where represents a perfect fluid and T^^ represents a scalar field with e = ±1. For e = +1, is regular and 
it represents the quintessence model. For e = — 1, the scalar field has a negative dynamic energy and Tj}^ represents 
the phantom model. And, we assume that each of the two components TJ^ and T^^, conserves independently just as 
in most of the 4D dark energy models. 

Here, we should notice that the coordinate time t is not the proper time in the solutions and, generally, one 
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cannot transform it to the proper time without changing the form of the 5D metric. However, on a given hypersurface 
y = const, the proper time r relates the coordinate time t via dr = Bdt. So, in this 5D model, the Hubble parameter 
H, the deceleration parameter q, and the statefinder parameter ([T]) introduced in Ref. ([23]) should be given as 
follows: 

H = 



(7) 
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Using the relations ([5])-® and the conservation laws T^™'*'' = and T^fi'' = 0, we can obtain the 4D Einstein equation 
and the equations of motion for the scalar field: 

Prn + 3 — (p,„ = 0, 

= - 4 ^iP"^ + +P<P+ P<l>)- (8) 

where 

1 (/)2 1 (/,2 

P<P^i^e^ + V{cly), p^^-e^~V{^). (9) 

are the energy density and pressure of the scalar field respectively. The potential V{(j)) is assumed to be exponentially 
dependent on by V{4>) — Vq exp(— Akc/)). pm and Pm are the energy density and pressure of the perfect fiuid 
respectively, and Pm = (7™ - l)pm.- 

The equation of state parameter for the scalar field is found to be: 

= ^ = £t^£!m. (10) 

III. STATEFINDER PARAMETER FOR 5D ATTRACTOR SOLUTION OF QUINTESSENCE MODEL 



The quintessence model is the scenario of equations (j8|)- (jl0|) with e = +1, which has been discussed in previous 
work [231. As in Ref. [2^, we define the x and y in a plane-autonomous system as: 

y = ^sv7- (11) 



For a spatially flat universe (fc = 0), we find that the evolution equation for x and y are of the same form as in Ref. 

/3" 3 
-y2 + -x[2x^ + 7„,(1 - x2 - y^)], 

y = -xJ\xy+^y[2x' + -ira{l-x^ -y"")], (12) 
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FIG. 1: The left figure is s — r diagram of attractor solution in quintessence model. The right figure is g — r diagram of attractor 
solution in quintessence model. The curves evolve in the variable interval N € [—2, 10]. Selected curves for jm = 1, A = 1, 1.2 
and 1.5 respectively. Dots locate the current values of the statefinder parameters. 



where a prime denotes a derivative with respect to the logarithm of the scale factor, N = hiA. And we can get the 
densities of two components: 

2 2 
K pjyi i22 r> P4> 2 , 2 /'10^ 

^ni = = ^-x^ -y , ^4'=' =^ +y ■ (13) 

and the equation of state parameter for quintessence is: 

? 2 

^'t' = 2 I 2 - (14) 

The statefinder parameters and deceleration parameter for this quintessence model are: 

r = (^7™- ^7m + l)(l-a;'-y')-3V6Axy2 + lOx2 + y^ 

^ (§7^ - f 7m + 1)(1 -x^- y^) - iV&Xxy'' + lOx' + y^ - I 
3[(|7™-l)(l-x2-y2) + 2a;2-y2„ 1] 

q = {\lrn-l){l-x^ -y'') + 2x^ -y\ (15) 

There are two late-time attractor solutions for which have been found in Ref. [23]. One is the scaling solution 
which has been discussed by the authors in Ref. 23], another is the late-time attractor solution dominated by the 
scalar field: 



X = A/%/6, y = Vl - AV6, with < 6 (16) 

In the following we will discuss the statefinder parameters for the attractor solution of quintessence model. In Fig. 
[T]we plot the evolution of the statefinder pairs {r, s} and {r,q} for the attractor solution of quintessence. The plot 
is for variable interval N € [—2, 10], and the selected evolution trajectories of r{s) and r(g) correspond to 7^ = 1, 
A = 1, 1.2 and 1.5 respectively. From the evolution trajectory of r(s) we can see that the curves pass through the 
fixed point LCDM in the past, while the distance from the curves to LCDM scenario is somewhat far in the future. 
The statefinder pair {r, s} lies in the regions s > 0, r < 1, which is differ from other quintessence model. The diagram 
of —ln[\ + z) — r and ~ln{l + z) — s in Fig. [5] show that the big deviation between statefinder parameter and 
the LCDM scenario has been caused with the increase of the parameter A in the potential V{(t>). From the diagram 
—ln(l + z) — g we can see that the deceleration parameter q satisfies g < at present. The Fig. [3]is the curves of the 
equation of state versus redshift with different A, which indicates that the universe is accelerating without appearance 
of the big rip in the future. 



IV. STATEFINDER PARAMETER FOR 5D ATTRACTOR SOLUTION OF PHANTOM MODEL 



The phantom model is the scenario of equations ([5)l- (|10p with e = — 1, which has been discussed in previous work 
1271. As in Ref. [271, we define the dimcnsionless variables: 
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FIG. 2: The figures are the statefinder parameters and deceleration parameter versus redshift diagram of the attractor solution 
in quintessence model. The left figure is — ln(l + z) — r diagram, the middle figure is — ln(l + z) — s diagram and the right 
figure is — ln(l + z) — q diagram. The curves evolve in the variable interval A'^ £ [—2, 10]. Selected curves for 7™ = 1, A = 1, 
1.2 and 1.5 respectively. Dots locate the current values of the statefinder parameters. 
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FIG. 3: The figures are the equation of state of the quintessence model w versus redshift. The curves evolve in the variable 
interval TV G [—2,5]. Selected curves for 7™ = 1, A = 1, 1.2 and 1.5 respectively. 

For a spatially flat universe (fc = 0), we find that the evolution equation for x and y in the phantom model are: 

x' = ^x[-2x^ +jrail + x'' -y^)]-3x- X^y^, 

y = ^y[-2x''+j^{l + x^-y^)]-\^xy. (18) 

where a prime denotes a derivative with respect to the logarithm of the scale factor, N = \nA. And we can get the 
densities of two components: 

n„.^^ = l + x'-y\ n,^^ = ~x' + y'. (19) 

and the equation of state parameter for phantom is: 

w, = (20) 
x'^ — 

The statefinder parameters and deceleration parameter for this phantom model are: 

r = - \l.n + 1)(1 + x'- y2) - 3V6Axy2 - lOx^ + y\ 

^ ihl^ - f7m + l){l + x^- y^) - SVeXxy^ - lOx^ + y' ~ 1 
3[(|7™-l)(l + 2;2-y2)_2a;2-y2_ 1] 

q = i^lrn - 1)(1 + x2 - y') - _ ^2^ (21) 

There is only one meaningful late-time attractor solution for (|18p . which is the scalar field dominated solution: 

X = -X/V6, y = V'TTAVe. (22) 
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FIG. 4: The left figure is s — r diagram of attractor solution in phantom model. The right figure is g — r diagram of attractor 
solution in phantom model. The curves evolve in the variable interval A'' G [—2, 10]. Selected curves for 7™ = 1, A = 1, 1.2 and 
1.5 respectively. Dots locate the current values of the statefinder parameters. 
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FIG. 5: The figures are the statefinder parameters and deceleration parameter versus redshift diagram of the attractor solution 
in phantom model. The left figure is — ln(l + z) — r diagram, the middle figure is — ln(l + z) — s diagram and the right figure 
is — ln(l + z) — q diagram. The curves evolve in the variable interval TV G [—2, 10]. Selected curves for 7™ = 1, A = 1, 1.2 and 
1.5 respectively. Dots locate the current values of the statefinder parameters. 



and the eigenvalue of this solution is (—3 — A^/2, —37m — A^). 

In Fig. m we show the time evolution of the statefinder pairs {r, s} and {r, q} for the attractor solution of the 
phantom field. The plot is also for variable interval N e [—2, 10], and the corresponding parameters are the same in 
the Fig. [TJ We also can see that the curves far from the LCDM in the future, but the statefinder pair {r, s} lies in 
the region s < 0, r > 1 which is different from that in Fig. [1] We see that the statefinder trajectories r — s is almost 
linear in the past and future, which means that the deceleration parameter changes from one constant to another 
nearly with the increasing of the time, which can also be seen from the Fig. [SI In Fig. O the curves of —ln{l + z) — 
r and —ln{l + z) — s show clearly that the parameter r will apart from 1 and s will far from with the increase of 
the parameter A. From —ln(l + z) — q diagram we can found that the deceleration parameter q of phantom are all 
satisfy q < —1 with different parameter A at present. In Fig. [Slwe see that the equation of state of phantom below 
— 1, which is different from the Fig. [31 and which could cause a big rip [ssl. [s^] in the universe in the future. 
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FIG. 6: The figures are the equation of state of the phantom model w versus redshift. The curves evolve in the variable interval 
N £ [—2,5]. Selected curves for 7™ = 1, A = 1, 1.2 and 1.5 respectively. 
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We have applied a statefinder analysis to the 5D attractor solutions of quintessence and phantom field, and investigate 
the effect of the parameter A in the potential V {(/)). Though these two attractor solutions are both scalar field 
dominated solution, the difference of these two solutions can be found in the statefinder parameters plane Fig[l]-Figl6l 
(i) The region of the statefinder pair {r, s} is different between Fig. [T] and Fig. U) For quintessence scenario Fig. [1] 
s > and r < 1 while for phantom scenario Fig. |4l s < and r > 1. (ii) The influence of the parameter A in potential 
V{(j)) is different in the evolution of the universe, which can be seen in the statefinder trajectories. In quintessence 
scenario Fig. [21 the increasing of the A will decrease the statefinder parameter r and the deceleration parameter q, 
while it will increase the statefinder parameter s. In phantom scenario Fig. 31 the statefinder parameter r will be 
increased, while the s and q will be decreased with the increasing of the A in the potential. So, through the statefinder 
parameter, we can see the difference between quintessence field and phantom field clearly. The attractor solution of 
quintessence field will cause the universe to accelerate forever, but, the attractor solution of phantom field will cause 
the universe to end with a big rip. So, We think that we should pay attention to find a suitable interaction between 
phantom and dark matter or other components in 5D model to avoid the big rip in our future work. 



V. DISCUSSION 



In summary, we have studied the statefinder parameter r and s in the five-dimensional big bounce model in a 
plane-autonomous system, and apply it to contrast the attractor solution of the quintessence field and the phantom 
field respectively. It is found that the evolving trajectories of these two attractor solutions in the {r, s} and {r, q} 
plane is quite different, and which is also different from the statefinder diagnostic of other dark models. Through the 
figures of the statefinder parameters and deceleration parameter we can see that the universe will accelerate forever 
in the quintessence scenario, while the universe will end with a big rip in the phantom scenario. We hope that the 
future high precision observation will be able to determine these statefinder parameters and consequently shed light 
on the nature of dark energy. 
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